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.,—1 . Abstract. Two-photon detachment from the F^ negative ion has been investigated 

^.^, within the lowest order perturbation theory. We show that in accordance with the 

T^ ' adiabatic theory a proper asymptotic behaviour of the 2p bound state wave function 

is crucial for obtaining correct absolute values of the multiphoton detachment cross 
sections. We find that the latter are substantially higher than it was previously 
beheved. 
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In the recent papers by Gribakin and Kuchiev (1997a, 1997b) an adiabatic analytical 
theory of multiphoton detachment from negative ions has been developed, based on the 
Keldysh approach (Keldysh 1964). Simple analytical expressions obtained there for the 
differential and total n-photon detachment cross sections allow one to estimate them for 
any negative ion. One of the important points of that work is that the electron escape 
from the atomic system in a low-frequency laser field takes place at large distances. 



r ~ 1/v/tJ ~ V2n//t > 1, (1) 

where u is the photon frequency, k is related to the initial bound state energy, Eq = 
— K^/2, and n is the number of quanta absorbed (atomic units are used throughout). 
Accordingly, the multiphoton detachment rates are determined by the long-range 
asymptotic behaviour of the bound-state wave function, namely by the asymptotic 
parameters A and k, of the bound state radial wave function R{r) ~ Ar~^e~'^^. This 
result has been obtained using the length form of the interaction with the laser field, 
which proves to be the most convenient for multiphoton processes. 

The analytical adiabatic approach is valid for multiphoton detachment processes, 
i.e., strictly speaking, for n ^ 1. However, the calculations for H~ and halogen negative 
ions indicate (Gribakin and Kuchiev 1997a, Kuchiev and Ostrovsky 1998) that the 
analytical formulae should give reasonable answers even for n = 2. The aim of the 
present work is to verify these conclusions by doing direct numerical calculations of 
the 2-photon detachment cross sections. In particular we examine the sensitivity of 
the photodetachment cross sections to the asymptotic behaviour of the ground state 
wavefunction and show that it is indeed very strong. Thus, a "small" 20% error in n 
present in the Hartree-Fock (HF) wave function of the fluorine negative ion results in a 
factor of three underestimation of the 2-photon cross section. This emphasizes the need 
to use bound-state wavefunctions with correct asymptotic behaviour in calculations of 
multiphoton processes. 

In this paper we calculate the 2-photon detachment amplitudes, cross sections 
and photoelectron angular distribution using the lowest-order perturbation theory and 
compare the results obtained with different ground state wavefunctions. We present and 
analyze the results for the F~ negative ion, where the results of a few other theoretical 
calculations (Robinson and Geltman 1967, Crance 1987a,b, Pan et al 1990, Pan and 
Starace 1991, van der Hart 1996) as well as experimental data (Kwon et al 1989, 
Blondel et al 1992, Blondel and Delsart 1993) are known. Pan et al and Pan and 
Starace calculated the two-photon detachment cross section and photoelectron angular 
distribution in the HF approximation (similar to that used by Crance) and with the 
account of first-order electron correlation effects. Their results show that the correlation 
corrections are about 15% for the partial and total cross sections and almost negligible 
for the angular distribution parameters, when the dipole length form is used. 



The total cross section of the two-photon detachment of an electron from an atomic 
system by a linearly polarized light of frequency oo is 
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where a^i is the partial cross section for the detachment into the final state with 
the photoelectron orbital momentum If and the total orbital momentum L, and the 
continuous spectrum wave function of the photoelectron is normalized to the 5-function 
of energy. For the 2p electron detachment from F^ 2p^ ^S the final state can be either 
^S {L = 0, If = 1) or ^D (L = 2, If = 1, 3). The two-photon amplitude Ai^l{uj) is 
determined by the following equations 
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where z// is the intermediate electron state with the orbital momentum / after absorbing 
the first photon (/ = 0, 2 for F~), and hqIq is the initial bound state. The reduced dipole 
matrix elements are defined in the usual way, e.g. in the length form 



(z//||ci||no/o) 
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where /> = max{/, Iq} and P's are the radial wave functions. 

The photoelectron angular distribution is described by the differential cross section 



da 

dn 



a 



An 



E/^2iMP2i(cOSI 



(6) 



3=0 



where 6 is measured with respect to the light polarization axis, and the asymmetry 
parameters f32j are defined in terms of the two-photon transition amplitudes Ai l and 
scattering phases of the photoelectron 5;^ : 
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where [/] = 2/ + 1 and /3o = 1, so that the photoelectron angular distribution after a 
two-photon detachment is characterized by P2 and P4. 

The self-consistent HF calculation of the F^ ground state yields the 2p-electron 
energy E^^ = —0.362 Ryd, which is much lower than its true value equal to the 



negative of the experimental electron affinity of F: E2p^ = —0.250 Ryd (Radtzig and 
Smirnov 1986). It is often assumed that the HF radial wavefunction is still a good 
starting point for calculations of multiphoton detachment, if the experimental binding 
energy is used in lieu of the HF value (Crance 1987a,b, Pan et al 1990). Pan et al 
(1990) showed that the two-photon detachment cross sections obtained with the dipole 
operator in the velocity form are very sensitive to the 2p-electron energy, while the 
length form results change little when the HF energy is replaced by the experimental 
one. However, one should use not only use the correct energy but, which is much 
more important, the bound state wavefunction with the correct asymptotic behaviourf. 
The importance of large distances, where one can use the correct asymptotic form of 
the bound-state wave function, speaks strongly in favour of using the length form of 
the dipole operator (Gribakin and Kuchiev 1997a). To correct the 2p wavefunction we 
solved the HF equations for the F~ ground state with an additional small repulsive 
potential V{r) = a/[2(r^ + a^)^]. We chose a = 1 and a = 0.61 au to ensure that 
the 2p energy was equal to the experimental value. The HF and corrected 2p radial 
wavefunctions P{r) = R{r)/r are presented in figure |l|(a). The difference between them 
appears to be small - it does not exceed 10% near the maximum. Their asymptotic 
behaviour P{r) ^ Aexp{—K,r) corresponds to A = 0.94 and 0.86, and k = 0.6 and 0.5, 
respectively. The difference in k means that the two wavefunctions are in fact quite 
different at large distances. 

The wavefunctions of the intermediate (ul) and final (sflf) states of the 
photoelectron are calculated in the HF field of the frozen neutral 2p^ core. The 
photoelectron is coupled to the core to form the total spin S = and angular momentum 
L: L = 1 for the intermediate / = 0, 2 states, L = 0, 2 for If = 1, and L = 2 for // = 3 
final states. The intermediate states continua are discretized and represented by a 70- 
state momentum grid with constant spacing Ap. 

There are two ways of calculating the two-photon amplitudes Mi^.i of equation @. 
The ffist is by direct summation over the intermediate states. It involves a non-trivial 
evaluation of the free-free dipole matrix elements together with the accurate treatment 
of pole- and S-tjpe singularities (Korol 1994, 1997). Another way of calculating such 
sums is by solving an inhomogeneous Schrodinger-type equation for the effective radial 
function of the intermediate state 

p ( s _X- Pui{r){ul\\d\\nQlo) 

^^'-^Eo + u-E^ + iO ^^^ 

(Sternheimer 1951, Dalgarno and Lewis 1955). This wavefunction describes the 

f The need for an asymptotically correct wavefunction is clearly illustrated by the adiabatic 
hyperspherical calculation of multiphoton detachment from H~ by Liu et al (1992), where a 3.4% 
change of k results in a 25% change of the 2-photon cross section. 



amplitude of finding the electron at different distances from the atom after absorption of 
the first quantum. After calculation of Puj{r) the amplitude is obtained from the radial 
integral as 

Mi^,{uj) = {-iy>^Jp,^,^.{r)P^{r)rdr , (/> = max{/^, /}) , (9) 

In the present work we calculate the two-photon amplitudes using both techniques. 
The second one is especially simple below the single-photon detachment threshold, 
uj < |ii^o|) where Puj{r) drops exponentially at large distances: Puj{r) oc exp(— i^^r), 
where K^ = [2(|i?o| ~^)Y^'^- At finite distances r < 1/Ap it can be computed easily by 
direct numerical summation in equation (|^). 

It is instructive to look at the shape of the effective wavefunction Pco{r). As an 
example, figure |l|(b) shows this function calculated for the intermediate d electron at 
uj = 0.226 Ryd for the HF ground state. The maximum of Puj{r) is shifted towards large 
radii, compared to the maximum of the ground state wavefunction. Also shown in figure 
|l](b) are the radial wavefunction of the final p-wave electron {sf = 2uj + E^^ = 0.09 
Ryd) and the integrand Pe^if{r)Pi^{r)r of equation @. These plots illustrate the point 
that the two-photon amplitude Mi i^^ (uj) is indeed determined by large electron-atom 
separations (||). Accordingly, the correct asymptotic behaviour of the ground state 
wavefunction is crucial. 

In this work we compare the cross sections and angular asymmetry parameters 
calculated in different approximations with both the HF and corrected 2p wavefunctions 
(figure |I|(a)). Let us first discuss the results obtained with the HF energy of the 2p 
electron. It corresponds to the two-photon threshold u = 0.181 Ryd. The cross section 
calculated from equations (0)-(|^) using the HF functions of the initial, intermediate 
and final states are shown in figure |^ by a short-dashed line. It is very similar to the 
dipole length lowest-order HF results of Pan et al (1990), although the latter are about 
10% lower than ours. What is the source of this discrepancy? Pan et al used the 
Roothaan-HF expansion of the bound state. This form of the bound-state wavefunction 
has an incorrect asymptotic behaviour at r > 7, see figure |l](a), inset. Because of the 
importance of large electron-atom separations in the multiphoton processes even a small 
error in the wave function could lead to some inaccuracies in the 2-photon detachment 
amplitudes. In the work of Pan et al electron correlation effects were calculated. It was 
shown that they suppress the cross section in F~ by about 20% at the maximum. These 
results shown in figure ^ by open squares are still close to the HF curve. 

The asymptotic parameters of the HF 2p wave function are k = 0.6016 and A = 0.94. 
We use them in the adiabatic theory formula (equation (5) of Gribakin and Kuchiev 
1997b) and obtain the cross section shown in figure Q by solid dots. It reproduces the 
energy dependence of the HF cross section well, though overestimates its magnitude by 
a factor of two. This is a reasonable result, since the adiabatic theory should only be 
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valid for n ^ 1. The calculations of Gribakin and Kuchiev (1997a) showed that for 
H^ and n = 3 the analytical adiabatic results are already 20% accurate. There are two 
approximations made in the adiabatic theory: (i) the use of the Volkov wavefunction 
to describe the photoelectron, and (ii) the saddle-point calculation of the integral over 
time, which enables one to express the amplitude in terms of the asymptotic parameters 
of the bound state. In the weak-field regime the use of the Volkov function is equivalent 
to the the so-called 'free-electron' approximation (examined earlier by Crance). In this 
approximation the wavefunctions of the photoelectron in the intermediate and final 
states are described by plane waves. When we do such calculation for F~ (chain curve 
in figure §) the results turn out to be very close to those of the adiabatic theory. This 
means that the approximation (ii) of the adiabatic theory is in fact quite good even at 
n = 2. 

When we use the experimental energy of the 2p electron together with the HF 
wavefunctions the magnitude of the 2-photon cross section changes very little (dashed 
line in figure ^, as seen earlier by Pan et al (1990) for both HF and correlated results 
(open circles). The HF results of Crance (1987a) are close to the above, and the cross 
section of van der Hart (1996) is also similar, with a maximum of 1.27 au at a; = 0.166 
Ryd. 

However, when we use the corrected 2p wavefunction, the photodetachment cross 
section increases more than three times. It is shown by solid line in figure ^ and 
we consider this to be the best evaluation of the cross section for F~. The cusp 
on the curve corresponds to the single-photon thresholdf. The same increase is also 
demonstrated by the adiabatic theory (with modified asymptotic parameters k = 0.4998 
and A = 0.86) and the plane-wave results. As we explained earlier, this "surprising" 
sensitivity of the multiphoton detachment probabilities to the asymptotic form of the 
bound-state wavefunction is a direct consequence of the dominant role of large electron- 
atom separations in this problem (Gribakin and Kuchiev 1997a). 

The error induced by the use of an asymptotically incorrect wave function can be 
estimated within the adiabatic approach. It turns out that if one uses the experimental 
binding energy Eq = —h?/2 together with an incorrect bound state P{r) oc exp(— kV), 
the n-photon cross section acquires an error factor 

u-.m\\ (10) 



where Ak = k' — n. This equation implies that the relative error in the amplitude is 
~ AkR, where R = I/a/cJ is the large radius from equation (|I]). For k' > n the error 
factor is smaller than unity. Thus, using a stronger bound wave function leads to an 



f This feature is a consequence of the Wigner threshold dependence a ex \^uj — Eq of the s-wave 
single-photon detachment from F^. 



underestimation of the cross section. For F~ the factor (|TOD calculated for Ak = 0.1 
and u = 0.085 au near the cross section maximum, yields 0.33. This value agrees with 
the difference between the cross sections observed in figure 0. The only other work that 
used an asymptotically correct 2p wave function was the model potential calculation of 
Robinson and Geltman (1967), which produced a cross section two times greater than 
those of Crance, Pan et al and van der Hart. 

To estimate the size of possible errors introduced by our way of correcting the 2p 
wavefunction we have examined the dependence of our cross section on the choice of a 
and a in the repulsive potential. We find that as long as the asymptotic behaviour of 
the 2p state remains correct, the two-photon cross sections are always proportional to 
that shown by the solid line in figure ^ Different pairs of a and a result in the variation 
of A, and the magnitude of the cross section is simply proportional to A^ . Our value of 
A = 0.86 is close to A = 0.84 from Radtsig and Smirnov (1986) and we are sure that 
our results are basically correct. Even a large 10% uncertainty in the value of A would 
mean a maximal 20% error in the cross section. In any case the cross section will be 
much larger than those obtained with the HP 2p ground state. 

The difference between experimental and HP values of the 2p energy is a 
manifestation of electron correlations. It influences the result via the asymptotic 
behaviour of the ground-state wavefunction. This is by far the most important 
correlation effect in multiphoton detachment. The use of the asymptotically correct 
2p wavefunction changes the cross section by a factor of three, which is much greater 
than other correlation effects (Pan et al 1990). This fact distinguishes this problem 
from the single-photon processes, where other correlation effects are essential. 

The angular asymmetry parameters (32 and /?4 calculated using the experimental 2p 
energy are shown in figure |^, together with the correlated length results of Pan and 
Starace (1991) and experimental points of Blondel and Delsart (1993) at u = 0.171 
Ryd. The asymmetry parameters (7) are relative quantities, and the results of 
different calculations are much closer for them than for the absolute values of the 
photodetachment cross sections. The adiabatic theory is again in good agreement with 
the plane-wave approximation, especially in f]^. It appears that this parameter is on the 
whole less sensitive to the details of the calculation, because it is simply proportional to 
the amplitude of / wave emission, and there is no interference in the sum in equation 
(7) for (54. The experimental values of f3 for P~ obtained in the earlier work of Blondel 
et al (1992) are close to those of Blondel and Delsart (1993). This is why P~ serves as 
a good benchmark for angular asymmetry calculations. The perfect agreement between 
adiabatic theory and the experiment is probably fortuitous. Pigure ^ indicates that for 
the experiment to be able to distinguish between various theoretical data one would 
wish to make measurements at higher photon energies, where the results of different 
approximations diverge. 



From the theoretical point of view it seems that the total cross sections and the 
angular asymmetry parameters are determined by different physical features of the 
problem. The absolute size of the cross sections is very sensitive to the asymptotic 
behaviour of the bound state wave function. This sensitivity increases for large-n 
processes, when uj become smaller, as suggested by estimate (0). The cross sections 
also depend on the atomic potential which acts on the photoelectron, hence the 
difference between the results obtained with the HF and plane waves. This latter 
effect is suppressed for larger n, since this re-scattering of the photoelectron is inversely 
proportional to some power of large R. The angular asymmetry parameters are affected 
by the electron-atom potential via the scattering phaseshifts 5^,. However, for large 
n and small photoelectron energy E r^ uj the phaseshifts should be close to integer 
multiples of vr. Besides that, contributions of higher partial waves become dominant. 
They are almost unaffected by the atomic potential and the adiabatic theory should 
become very accurate. 

In summary, we have shown by direct numerical calculations that in agreement 
with the adiabatic theory, the multiphoton detachment rates are very sensitive to the 
asymptotic behaviour of the bound state wavefunction. For fluorine this means that the 
true 2-photon detachment cross sections are substantially higher than it was believed 
earlier. The discrepancy revealed is much greater than other electron correlation effects. 
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Figure captions 



Figure 1. Wavefunctions of the F~ ground state, effective intermediate state and 
final state of the photoelectron. (a) Radial wavefunction of the 2p subshell of F~ in 

the HF approximation ( , fi'Sf = —0.362 Ryd), and that with a model potential 

added to reproduce the experimental energy ( — • — , E2p — —0.250 Ryd). The 
inset shows the same wavefunctions on the logarithmic scale, together with the radial 

Roothaan-HF 2p radial wavefunction of F^ from Clementi and Roetti (1974), . 

(b) , HF 2p wavefunction; , effective wavefunction P,^{r)^ equation (0), 

of the intermediate d state at a; = 0.226 Ryd; — • — , final state p Q^D) wavefunction, 
e = 0.09 Ryd; , integrand of equation (^ for the 2-photon amplitude Mpd- 



Figure 2. 2-photon detachment cross sections. Present calculations: , HF 

wavefunctions of the 2p, intermediate and final states; , same with the corrected 2p 

wavefunction; — • — , using plane waves in the intermediate and final state; , 

HF wave functions combined with the experimental 2p energy; •, adiabatic theory 
(equation (5) of Gribakin and Kuchiev 1997b). Other results: □ and o, calculations 
by Pan et al (1990) with the HF and experimental binding energies, respectively; ■, 
experiment (Kwon et al 1989). 



Figure 3. Photoelectron angular distribution parameters. Present calculation: 

, HF wavefunctions of the initial, intermediate and final states, experimental 

2p energy; , corrected 2p wavefunction, HF intermediate and final states; — • — , 

same with the plane wave in the intermediate and final states; •, /? parameters 
obtained from the adiabatic theory (equations (3), (4) of Gribakin and Kuchiev 1997b). 
Other results: o , correlated length results by Pan and Starace (1991); ■, experiment 
(Blondel and Delsart 1993). 
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